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Abstract 

We present a method, which constructs a suitable deformation vector in loop momentum 
space, when the loop integration is done numerically with the help of the subtraction method. 
The method presented here extends previously discussed techniques from the massless case 
to the general case of arbitrary masses in the loop. 



1 Introduction 



NLO predictions for multi-jet final states are important for the LHC experiments. In recent years 
there has been tremendous progress in our ability to compute observables at NLO with a sizable 
number of final state partons. NLO calculations with six or seven final state partons are now 
feasible [IIH31. One of the techniques which allow to go to such high parton multiplicities is a 
purely numerical approach. Within this approach also the virtual one-loop corrections are treated 
by Monte Carlo methods. In the numerical approach the NLO contribution to an observable 
is given as the sum of three sub-contributions. The first one is the real emission part minus 
suitable and well-known subtraction terms ll4]-[T0ll. The subtraction terms ensure that the phase 
space integral can be performed in four space-time dimensions with Monte Carlo techniques. 
The computation of this part is similar for all groups working in the field of multi-parton NLO 
calculations. 

The techniques of the various groups differ how they compute the virtual part. Within the 
numerical approach the ideas of the subtraction method and Monte Carlo integration are carried 
over to the virtual part. One subtracts from the one-loop amplitude suitable approximation terms 
for the soft, coUinear and ultraviolet singularities [fTTUTSl . The difference is integrable and the 
integration over the loop momentum can be combined with the integration over the phase space 
of the final state particles in one Monte Carlo integration. 

The third contribution to an NLO observable is given by the integrated subtraction terms from 
the real and from the virtual part. The sum of the real and virtual subtraction terms is finite and the 
phase space integration can again be performed by Monte Carlo methods. From a computational 
point of view this contribution is only slightly more expensive than the corresponding Bom 
contribution. 

The essential part of the numerical method is the treatment of the virtual corrections. Apart 
from the subtraction terms for the integrand of the one-loop amplitude the second important 
ingredient is a method for the contour deformation. In the vicinity of a singularity we use, 
whenever possible, contour deformation to by-pass the singularity. In the case where this is not 
possible, because the contour is pinched, and if the singularity is not integrable then (and only 
then) there is a subtraction term for it. The quality of the choice of the integration contour directly 
translates into the final Monte Carlo integration error. 

Several algorithmic approaches for the contour deformation have been discussed in the liter- 
ature [[T2] - l22l . They can be classified roughly according to three categories: In the first category 
one introduces Feynman parameters and integrates (after Feynman integration) over the loop mo- 
mentum analytically. Contour deformation (possibly in combination with sector decomposition) 
is applied to the integration over the Feynman parameters. Also in the second category Feynman 
parameters are introduced, however the loop integration is not done analytically but numerically. 
Contour deformation is applied to the integration over the loop momentum and the integration 
over the Feynman parameters. Within the third category no Feynman parameters are introduced 
and contour deformation is applied to the integration over the loop momentum. This is called 
"direct contour deformation". An algorithm to construct the contour deformation within the third 
category has been given in ref. lfT6]| for the case where all internal particles are massless. 

In our calculation for — 7 jets it turned out that the direct deformation method performs 
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best. In view of the LHC experiments, where massive particles like the Higgs, the top quark, the 
W- and Z-boson or hypothetical heavy particles beyond the Standard Model play an important 
role, it is desirable to generalise the direct deformation method to arbitrary masses. This is the 
purpose of the present paper. We present an algorithm, which constructs the deformation vector 
in the general case. In particular we can have some internal masses non-zero and others equal 
to zero. We assume that if an internal mass is non-zero then it is non-negligible against the 
centre-of-mass energy of the experiment. This is the case for heavy particles like the Higgs, 
top, W- and Z-bosons at the LHC. It excludes the treatment of mass effects of extremely light 
particles like the electron in LHC experiments (which also for other reasons are non-trivial to 
treat numerically). 

The generalisation to arbitrary masses is not a straightforward task. The method of ref. lfT6ll 
for the massless case is built on the fact that the singularities are located on light-cones, which are 
nested inside each other. Unfortunately the property that the singular surfaces are nested inside 
each other no longer holds if we replace the cones by more general hyperboloids. We solve this 
problem with two techniques: On the one hand we introduce additional reference points within 
the intersection of forward with backward mass hyperboloids. On the other hand we probe pre- 
defined standard directions in order to check if a deformation in this direction is possible. 

This paper is organised as follows: In section [21 we define the problem and introduce our 
notation. Section[3]is the main part of this paper and specifies the algorithm for the contour inte- 
gration. In section m we describe checks which we have performed to test our method. Finally, 
section [5] contains our conclusions. 



2 Numerical integration of loop amplitudes 

In this paper we consider primitive one-loop amplitudes. Primitive amplitudes are gauge-invariant 
building blocks of the full one-loop amplitude. Each primitive amplitude has a fixed cyclic order- 
ing of the external particles. The cyclic ordering ensures that there are at maximum n different 
loop propagators in the problem, where n is the number of external legs. We label the external 
momenta clockwise by p\, pi, Pn and define qi = pi+ P2 + ■■■+ Pi- The loop momenta are 

j 

1=1 

For convenience we set 

ko = kn and qQ = q„. (2) 

Due to momentum conservation we actually have 

<?o = <?« = 0. (3) 

Nevertheless we will use (or qn), this makes the formulae more symmetric with respect to the 
indices. The labelling of the momenta is shown in fig. ([T])- 
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Pi 



Pn 



Figure 1 : The labelling of the momenta for a primitive one-loop amplitude. The arrows denote 
the momentum flow. 



In addition we define an ultraviolet propagator with momentum 

k = k-Q (4) 

and mass squared /j^y. The four-vector Q and the parameter jj^y specify the ultraviolet prop- 
agator. The mass squared /j^y can be chosen freely, and a useful choice is to take ^u^v Purely 
imaginary with Im ^u^y < 0. Choosing the value along the negative imaginary axis large enough 
ensures that the ultraviolet propagator never goes on-shell. 
The object of investigation is an integral of the form 

^ f d^k Rjk) 

where is a rational function of the loop momentum k^, which has only poles at P —f^y = 0- 
With a proper choice of /^uv ensure that the integration contour always stays away from 

the poles of the ultraviolet propagator. Therefore the only poles which we have to care about 
are the ones which are shown explicitly in eq. (O. It is assumed that the integral is finite, i.e. 
all contributions which would lead to poles in the dimensional parameter £ when the integral 
is calculated in D = 4 — 2£ dimensions are absent. This can be achieved with the subtraction 
method. Usually the integrand is the integrand of a primitive bare one-loop amplitude minus 
subtraction terms for the soft, coUinear and ultraviolet singularities. Suitable subtraction terms 
can be found in the literature lfTTl[T2l[T4l . Since the integral is finite, it can be computed in 
four dimensions. However, this does not yet imply that we can safely integrate each of the four 
components of the loop momentum k^ from minus infinity to plus infinity along the real axis. 
There is still the possibility that some of the loop propagators go on-shell for real values of 
the loop momentum. If the contour is not pinched this is harmless, as we may escape into the 



4 



complex plane in a direction indicated by Feynman's +/5-prescription. However, it implies that 
the integration should be done over a region of real dimension 4 in the complex space C^. If the 
contour is pinched then the singularity is integrable when the integration is done over the loop 
momentum space and the phase space. This is the case because either the singularity in the bare 
one-loop amplitude is integrable by itself, or - if not - there is a subtraction term for it. 

The purpose of this paper is to give a method for a suitable deformation into the complex 
space C^. To this aim we set 



k = k + i}c(k) , 
where kf^ is real. After this deformation our integral equals 



/ 



d'^k 
(2^ 



dkf 



mm 



The Jacobian 



n (kj-mj-K^ + 2ikj-K 
7=1 ^ 



dk^' 



(6) 



(7) 



dk" 



(8) 



for the contour deformation can be computed numerically. To match Feynman's +/5-prescription 
we have to construct the deformation vector K such that 



kl 



m.; 







— > 



kj-K> 0, 



(9) 



and the equal sign applies only if the contour is pinched. As long as K is infinitesimal we may 
neglect the term in eq. ([7]). On the other hand in order to minimise Monte Carlo integration 
errors we would like to choose K as large as possible. We will therefore first construct the 
direction of K and then scale K to the maximal length such that no poles are crossed by varying K 
from zero length to its maximal length. 

The singularities of the integrand lie on mass hyperboloids defined by {k - 
origins given by qq, qi, ^„_i. Since 



— m? = with 



(10) 



the external momenta p, connect the origins of the mass hyperboloids and we arrive for a generic 
primitive one-loop amplitude at the graphical representation shown in fig. Q on the left. The 
lines connecting the points qi form a closed loop due to momentum conservation 



(11) 



i=\ 



Without loss of generality we may assume that one of the two momenta corresponding to the 
initial-state particles is p„. We then denote the momentum of the other initial-state particle by 
Pa- In the case where the two initial-state particles are adjacent we can take a = n—l and the 
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Figure 2: The origins of the mass hyperboloids in loop momentum space. The diagram on the left 
shows the generic case for a primitive amplitude, the diagram on the right shows the degenerate 
case, where the two incoming particles are adjacent. 



diagram degenerates to the one shown in fig. ^ on the right. In practice the two initial-state 
particles are always massless and the two initial-state momenta are light-like: Pn = Pa = ^■ 
We define the forward mass hyperboloid as all points, which satisfy 

{k-qif-mj = 0, k^>ql (12) 
Analogously, the backward mass hyperboloid is defined by 

{k-qif-mj = 0, k^<ql (13) 

The interior region of a mass hyperboloid is given by the set of all points, which satisfy 

{k-qif -m}>0. (14) 

We use the convention that the boundary is part of the interior region of a mass hyperboloid. 

In loop momentum space the points ^i^ •••j Qn-i are contained within a finite region. We 
call this region the "interior region of the loop momentum space". If it is clear from the context 
we will just use the term "interior region" for the interior region of the loop momentum space. 
We may define the interior region of the loop momentum space by the set of all points, which 
are in the interior of some forward mass hyperboloid and in the interior of some backward mass 
hyperboloid. The exterior region of the loop momentum space is defined by the complement 
of the interior region of the loop momentum space. We note that only in the interior region 
intersections of forward hyperboloids with backward hyperboloids occur. 

3 Contour deformation 

In this section we show how to construct the deformation vector K. We write 

K = X(Kint + Kext), (15) 
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where Kint is a deformation vector for the interior region of the loop momentum space, Kext is a 
deformation parameter for the exterior region of the loop momentum space and X is a scaling pa- 
rameter. We discuss the construction of the two contributions Kext, Kint and the scaling parameter 
X separately. 

Some functions in the construction below will depend on individual components of four- 
vectors and are therefore dependent on the Lorentz frame. To make the definition unique we 
choose a specific Lorentz frame. A possible choice is the centre-of-mass system of the event. We 
make the convention that all frame-dependent statements refer to the centre-of-mass frame. 

In order to make this section more readable we will write throughout this section for the real 
loop momentum k instead of k. 

3.1 Helper functions 

We define some functions which are used in the construction of the deformation vector. We start 
with two functions h^_^ and defined by 

(±k° - VP + mA 

h^{k,m\Mj) = / . (16) 

(±kO-VP + mA +M\ 



The interpretation of the two functions /zg^ is as follows: The function /z§_|_ vanishes whenever k 
lies on the forward mass hyperboloid defined by 



a „.2 _ Q k^ > 0. (17) 



m 



Analogously the function vanishes whenever k lies on the backward mass hyperboloid de- 
fined by 

k^-rn^ = and k^ < 0. (18) 

If k is far away from the hyperboloid the function converges to /zg^ = 1 . The parameter 
controls the size of the transition region, where the function increases from zero to one. A 
typical value for the parameter is given by Mi = 0.035 where a/? is the centre-of-mass 
energy of the partonic event. 

We further define a function /z§ by 

(\k^\ - VP + mA 
h^{k,m\Ml) = J . (19) 

(\k^\-^JP + m^) +Mj 

The function /i§ vanishes whenever k lies on the (forward or backward) mass hyperboloid defined 
by 

k^-rn^ = 0, (20) 
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and converges to one far away from the hyperboloid. The functions h^j^_ and /?§_ are similar to 
the functions and /z defined in lfT6ll for the massless case, except that they do not contain 
any Heaviside step function. As a consequence they only cut out a region near the forward or 
backward mass hyperboloid, respectively, but not the region inside the mass hyperboloid. 
We further define a function /ze(f,Mj) by 

h{tM) = J^^^')^ (21) 

where 0(f) is the Heaviside step function. The function hQ{t,Ml) vanishes whenever the variable 
t is smaller zero and converges to one for f ^ 0. 

3.2 The deformation in the exterior region 

We start with the definition of the deformation vector Kext for the exterior region. In the exterior 
region we have per definition no intersections between forward and backward mass hyperboloids. 
Therefore we can separate the cases, where the loop momentum lies on a forward hyperboloid 
from the cases, where the loop momentum lies on a backward hyperboloid. Pinch singularities 
do not occur in the exterior region. As a further requirement on the deformation vector in the 
exterior region we have to demand that the ultraviolet power counting is respected. In other 
words we have to ensure that each propagator falls off like k^^ for k ^ oo independent of the 
direction in which we approach infinity. 

We write the deformation vector Kext as a sum of two terms. The first term deforms correctly 
if the loop momentum lies on a forward hyperboloid, while the second term deforms correctly 
if the loop momentum lies on a backward hyperboloid. The deformation vector for the exterior 
region is given by 

<xt(^) = g.vic+kl + C^k'L) ill) 

with 

k± = k-P±. (23) 
and g/jyj = diag( 1 , — 1 , — 1 , — 1 ) the metric tensor. The four- vectors P± are given by 

P+ = Z+{qa + qo,qa-qo), 

P = Z-{qa-i+qn-uqa-l-qn-l), (24) 

with 

Z'^{x,y) = ^-^±|UV-^V)). (25) 

The four- vector P+ lies on the intersection of the backward light-cone from with the backward 
light-cone from qa in the plane spanned by the four-vectors qo and qa- The four-vector P lies 
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<^4 



Figure 3: The definition of the four- vectors / + and P in the loop momentum space for a primi- 
tive amplitude with n = l and a = 4. 

on the intersection of the forward light-cone from qa-i with the forward light-cone from ^„ i in 
the plane spanned by the four- vectors qa-i and q,i-i- This is shown graphically in fig. ([3]). 
The coefficients c± are defined as 



i=l 



If the loop momentum k lies on a backward hyperboloid defined by 

kf-mj = and < 0, ie{l,...,n}, 



(26) 



(27) 



then one has h^_{ki,mj,M^) = and therefore c+ = 0. Analogously, if the loop momentum k 
lies on a forward hyperboloid given by 



m,- 







and 



k^ > 0, ie{l,...,n} 



(28) 



then we have /z§^(fc,,m^,Mj) = and therefore c_ = 0. The coefficients c± ensure that all 
contributions vanish which would contribute with the wrong sign to the imaginary part. It is easy 
to show that Kext always deforms correctly in the exterior region. Let us assume again that the 
loop momentum k lies on a forward hyperboloid as in eq. (l28l) . It follows that in this case 



k? > 



ki 



(29) 



At the same time k lies inside the forward cone with origin and we also have 



> 



(30) 
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We further have in this case C- = 



and c+ ~ 1 . For the imaginary part we have to consider 



Kext ■ k ^ k+o ki > 



(1+COS0+,/) >0. (31) 



Here we denoted by o the Euclidean scalar product of two four-vectors and by 0+,, the angle 
between the spatial three- vectors k+ and ki. 

We now discuss the ultraviolet behaviour of the deformation. We note that kf can be small 
even if the loop momenta k is in the ultraviolet region. We demand that the propagator falls off 
like I /{k ok) in the ultraviolet region after contour deformation. In the exterior region at least one 
of the coefficients c± is unequal to zero. The scaling parameter X is of order 1 in the ultraviolet 
region. The imaginary part of a propagator in the ultraviolet region can be estimated by 

ki-K ^ ^-Kext ^ kok, (32) 

and therefore all propagators fall off like l/{kok) as desired. 



3.3 The deformation in the interior region 

We now consider the definition of the deformation vector Kint for the interior region. Here we face 
several situations: In the interior region we have to take regions into account where a forward 
hyperboloid intersects with a backward hyperboloid. There can be regions which are bounded 
by more than two hyperboloids. If massless particles are present we can also have cones which 
are tangential to each other. We split the deformation vector for the interior region into three 
parts. The first part is sufficient for the massless case, while the second part provides a defor- 
mation when two non-zero mass hyperboloids intersect. Intersections of three or more mass 
hyperboloids are handled by the third part. We write 

<„t = I ^<7^ + <oft (33) 

i=i i,j=i 

><j 

with 



ki = k-qi and kij = k-Vij. (34) 
We note that if all particles are massless it is sufficient to define 

<nt = (35) 

(=1 

Up to minor modifications the definition in eq. (l35l) corresponds to the one given in ref. [fT6l . 
However, if massive particles are present the extra terms in eq. (I33l) are required. We start by 
giving the definitions of the coefficients c, and Cij, then discuss the vectors v/y and specify Kjoft in 
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the end. The coefficients are defined such that they are zero if the corresponding vector deforms 
into the wrong direction. In formulae 



m, 



and 



—Ciki -k] < 
-Cijkij -ki < 



Ci 



0, 
0, 



V/G (36) 



This ensures that the deformation vector Kint never deforms in the wrong direction. We set 



g (^centre , Yl , ^1) 11 ^^'^ 

1=1 



g(Wre,Yl,Mf)fJj/^-/' 
1=1 



(37) 



with the function 



YlM; 



kok + MV 



^centre — 2 



(38) 



The momenta k^ and k- have been defined in eq. (1231) . The function g ensures that the deforma- 
tion vector Kint fells off like l/{kok) in the ultraviolet region. Typical values for the parameters 
Yi and A/f are yi = 0.7 and M2 = 0.7 Vi. 
The factors dij are defined by 



d, 



i.l 



h8+{ki,mj,Ml 
h^_{ki,m],M\\ 
^ max [h^{ki,m],Ml)M{-^h-kuMl 



I = i, mi = 

{qi-qif = 0, q^<q^. m; = 
{qi-qif = 0, q^i>q^l, mi = 
otherwise 



(39) 



The factors djjj are defined by 

dij,i = hQ{zij,Mj) msix[h8{ki,mj,Mj),hQ{-2krkij,Mj 
Zij = 



[qi-qjf - {mi + nijY . 



(40) 



Let us explain the idea behind these definitions. The first three lines in eq. ( |39l ) handle exceptions 
due to massless particles in the loop. In coUinear configurations the contour is pinched in all 
directions perpendicular to the coUinear axis, but not along the coUinear axis. The first three 
lines allow for the possibility to deform along the coUinear axis. 

In all other cases one considers the conditions kf = mj and if —ki -ki < 0. If both condi- 
tions are fulfilled djj is zero. This ensures that all terms contribute with the correct sign to the 
imaginary part. 

The four-vectors Vij are defined by 



2 I 



(mj-mj) 
11 



(41) 




Figure 4: Definition of the vector v,j. The line connecting the points cjj and qj cuts the mass 
hyperboloids at the two points Vij and Vji. The vector Vij is at half of the distance between these 
two cutting points. 




m 



" k' 

Figure 5: If three or more mass hyperboloids intersect, the vector vij need not lie inside the 
intersection region. 
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The geometric interpretation of the vector v,j is shown in fig. dH). 

If three or more mass hyperboloids intersect, it can happen that none of the vectors v/y lies 
inside the intersection region. An example of this situation is shown in fig. ©. In this case the 
vectors defined up to now do not provide a deformation, since all coefficients vanish. The full 
interior region has a typical length scale of v^. Regions, where three or more mass hyperboloids 
intersect have typically a much smaller length scale. The vector Kgoft provides a deformation in 
these regions. Since the algebraic formulae describing the intersection of three or more mass 
hyperboloids are already rather complicated to handle, we use for the construction of Ksoft a 
different approach. In order to find a suitable deformation vector we probe a pre-defined set of 
directions and accept a given direction, if the deformation goes in the right direction. We write 

Ksoft = J^CaKa, (42) 

a 

where the index a sums over all pre-defined directions. It turns out that the set of the four 
directions along the Cartesian coordinate axis does already a good job. We set 

Ko = £soft (1,0,0,0), Ki = £soft (0, 1,0,0), 

K2=£soft(0,0,l,0), K3=£soft(0,0,0,l), (43) 

with £'soft being an energy scale much smaller then the centre-of-mass scale. A typical value is 
£^soft = 0.03v^. The coefficients c" are given by 

Ca = ^(We,Yl,M2)(^nJ+-nJ,-,j (44) 

with 

= max[/z5(fcz,m^,y2Mf),/ze(2fc,-K'',y2M?)] , 
d~i = max[/?5(/r/,m^,Y2M?),/ze(-2fc/-K^y2M?)]. (45) 

The product of the J^/s probes, if a deformation in the direction of is allowed, while the 
product of the d^j's probes if a deformation in the direction of (— k^) is allowed. Far away from 
all mass shells all functions d^^ are one and the contribution from Kq vanishes. If we approach 

a mass shell, the values of the functions d^i start to differ, thus yielding the right deformation 
vector. Note that the factor 72 in eq. (|45] ) ensures that also the transition region for the functions 
and has a smaller scale. A typical value for 72 is 72 = 0.008. 



3.4 The scaling parameter X 

Having defined the deformation vectors Kext and Kint for the exterior and interior region, respec- 
tively we now discuss the definition of the scaling parameter X entering the definition of the full 
deformation vector k: 

K = X (Kint + Kext). (46) 
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For the sum of Kjnt and Kext we write 



Ko 



Kint + K 



ext- 



(47) 



The definition of the scaling parameter X is along the lines of [16] with some minor changes due 
to the masses which can appear in the loop propagators. 

In the previous sections we discussed the construction of Kq = Kint + Kext in detail. By con- 
struction Kq deforms in the right direction. It follows immediately that if X is infinitesimal we 
have a correct deformation. But for small Monte Carlo errors we would like to deform as far as 
possible away from the singularities and therefore we need to make X as large as possible. In 
doing so we have to ensure that we do not cross any poles by varying the size of X from zero to 
its final value. 

We define for each propagator in the loop a scaling parameter Xj such that we do not cross 
any pole of the given propagator by varying Xj from zero to its final value. By taking A, as the 
minimum of these Xj we ensure that we do not cross any poles. We write for the j'-th propagator 



[kj + iXko)^ - mj = kj + 2iXKQ-kj-X^K^-mj 



The function Dj vanishes for values of X given by 

X 

where we introduced the functions 



Kn ■ ; / 

^0 



Ko ■ kj 



If Yj > Xj and Xj — we have a pole at 



and Yi 



Yj- 



(49) 



(50) 



(51) 



To avoid this pole we limit the value of X to one half of this value. If instead Yj < Xj the poles of 
the propagator Dj are on the imaginary X-axis and we can choose the real value of X as large as 
we like. We define 



Yj/4 
Xj-Yj/4 
Xj-Yj/2 



2Xj < Yj 

<Yj < 2Xj 



(52) 



Next we define a value Xcoii by 



with 



X, 



coll 



Yj < 



4C' 



(53) 



(54) 
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?icoii has two functions: First of all, in the case where several coefficients from the set of the c/ 
and Cij are non-zero, the scaling parameter Xcou acts as an averaging procedure. Secondly, in the 
collinear case the numerator and the denominator in the definition of Xj in eq. (l52l) goes to zero 
and A-coii ensures numerical stability in this limit. 
In addition we define a parameter ?iuv by 

f 1 : 4Ko-fc > lm{^ly) 

XuY ensures that the imaginary part of the ultraviolet propagator is always positive. 
The scaling parameter X is then given by 

X = min[l,Xi,...,X„,XcoiiAuv] ■ (56) 



4 Checks 

In order to test our method we verify that well-known one-loop integrals are computed correctly. 
We fix n external momenta pi, p2, Pn and n internal masses mi, m2, We let Sr be a 

subset of {!,...,«} containing r indices and we consider the integral 



/ 



,4, llikf-mj) 



(2^)^ n 



The integral is a scalar one-loop {n — r) -point function 

(271 



, d^k 1 

I - T^^^r^ r, (58) 



U{kj-mj 



and we consider scalar one-loop {n — r) -point functions which are infrared and ultraviolet fi- 
nite. We compare our results with the results obtained from the program "LoopTools" [l23ll . The 
contour deformation is obtained from the n propagators in the denominator of eq. (1571) and inde- 
pendent of the numerator in eq. (l57l) . Although the analytical result is given as a (n — r)-point 
function, the contour deformation corresponds to a n-point function. This is illustrated in fig. Q. 
By varying the set Sr in the numerator we can probe different regions in loop momentum space. 

We have verified the correctness of our method for randomly chosen external momenta and 
various choices of the internal masses. We have found good agreement with the known analytical 
results at the per cent to the per mille level for a typical Monte Carlo integration with 20 ■ 10^ 
evaluations. The most challenging parts from a numerical point of view are thresholds. For 
example, one-loop functions may contain terms of the form 

\n^{-s + m^). (59) 
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Figure 6: The diagram on the left defines the contour deformation for a process with n — 2 
massless and 2 massive external legs. Pinching {n — 3) loop propagators results in a three point 
function shown on the right. 
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Figure 7: Comparison of the results obtained by Monte Carlo integration with the analytical 
results in the vicinity of a threshold. 
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Terms like these will give at the threshold s = discontinuities and logarithmic divergences. 
In fig. (|7]) we show a typical situation, where the real part exhibits a discontinuity while the 
imaginary part shows a divergence. The analytical results have been obtained with "LoopTools" 
ll23l . while the numerical results have been obtained by Monte Carlo integration with the method 
described in this paper. The figure clearly shows, that also in the vicinity of thresholds the 
numerical method gives reasonable results. The deformation defined by Kgoft plays an essential 
part in the vicinity of a threshold. 

5 Conclusions 

The numerical calculation of the virtual corrections rely on a method for the contour deformation. 
In this paper we have given an algorithm to construct the deformation vector in loop momentum 
space for arbitrary masses in the loop. The algorithm presented in this paper opens the possibility 
to treat massive particles within the direct deformation method. 
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